Physisorption of positronium on quartz surfaces 
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The possibility of having positronium (Ps) physisorbed at a material surface is of great fundamen- 
tal interest, since it can lead to new insight regarding quantum sticking and is a necessary first step 
to try to obtain a Ps2 molecule on a material host. Some experiments in the past have produced 
evidence for physisorbed Ps on a quartz surface, but firm theoretical support for such a conclusion 
was lacking. We present a first-principles density-functional calculation of the key parameters deter- 
mining the interaction potential between Ps and an a-quartz surface. We show that there is indeed 
a bound state with an energy of 0.14 eV, a value which agrees very well with the experimental esti- 
mate of ~ 0.15 eV. Further, a brief energy analysis invoking the Langmuir-Hinshelwood mechanism 
for the reaction of physisorbed atoms shows that the formation and desorption of a Ps2 molecule in 
that picture is consistent with the above results. 

PACS numbers: 68.43.-h, 36.10.Dr, 34.50.Dy 



o 



> 

cn 
o 
p 

in 
o 
1^ 
o 



X 



Since positronium was first observed in 1951,^ positro- 
nium physics has evolved into a vigorous and fascinat- 
ing field, at the intersection of nuclear, atomic and con- 
densed matter physicsj2 The positronium (Ps) atom, i.e., 
an electron-positron bound state, is rather unique in that 
it is one of the very few examples of exotic particle- 
antiparticle atoms that have been observed, and the only 
one in the lepton family (the other examples, such as 
charmonium, being in the hadron family) . Ps is produced 
by irradiating matter with positrons, and its formation 
and stability has been studied in a variety of solid hosts.'' 
In the case of metals, the injected positrons thermalize 
and can reach the surface, where they can bind with an 
electron to form Ps. Thus, Ps desorbed from Cu and Al 
surfaces has been well measured and characterizedi^ In 
the case of insulators, in which Ps can also be formed 
in the bulk, Ps emission from AI2O3, Si02, and MgO 
powders was reported even earlier 

A question of great interest in these systems is whether 
Ps can be bound to the surface of the host. Indeed, this 
is related to fundamental questions such as the nature 
of quantum sticking^ and the potential of a Ps trap to 
lead to dipositronium (PS2) formation and its connection 
to the realization of a Ps Bose-Einstein condensates^!* In 
1976, intriguingly low decay rates in precision Ps lifetime 
measurements in Si02 powders were reported.- These 
findings were analyzed theoretically considering the pos- 
sibility of Ps physisorption on the surface of the grains, 
but the conclusion was that no surface bound state 
was possible. Later, however, temperature-dependent Ps 
emission measurements produced experimental evidence 
suggesting Ps physisorption on an a-quartz single crys- 
tal surface, with an estimated binding energy of ~ 0.15 
eV. Very recently, moreover, Ps lifetime measurements of 
high-density Ps gas in nanoporous silica suggest the for- 
mation of Ps2 molecules, a conclusion assuming that Ps 
atoms are previously physisorbed on the surface of the 
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Consequently, we think that it is important to re- 
consider the possibility of Ps physisorption on quartz 
surfaces. In the case of Al surfaces, for instance, the 
longer than expected (again) positron lifetimes observed 
in experiment^ can be understood in terms of a model 
calculatior^ii of weakly physisorbed Ps on the metal 
surface. Although an alternative explanation invoking 
trapped positrons in a slowly varying inhomogeneous 
electron gas was proposed by other authors, the for- 
mer possibility still exists.^*' Moreover, it is important 
to recognize that the previous theoretical estimate men- 
tioned was only approximate.^*' It is clear that a more 
accurate calculation of the Ps-quartz surface interaction 
potential is essential in order to discard or support the 
Ps physisorption interpretation of the above experiments 
and to find reference values for the physical parame- 
ters characterizing the Ps-surface interaction. Atomic 
physisorption on single-crystal surfaces has been exten- 
sively studied in the past. An account of the typical 
values found in experiment several characteristic param- 
eters can be found in Ref. [l3- From the theoretical point 
of view, following the pioneering work of Zaremba and 
Kohn,^* a good calculational scheme emerged thanks to 
the contribution of several authors, which we essen- 
tially follow here. In this respect, of foremost importance 
are accurate calculations of the dielectric function of the 
crystal and of the electronic charge density profile at its 
surface. In this work, we present a density-functional 
calculation of these two quantities. Our results for the 
Ps/a-quartz surface interaction potential show that in- 
deed Ps has surface bound states. As will be seen be- 
low, we find that the ground state binding energy is very 
close to the experimentally suggested value, which is also 
consistent with the possibility of of Ps2 formation via 
the Langmuir-Hinshelwood reaction mechanism for ph- 
ysisorbed atomsi^ 

It is well established that the interaction potential be- 
tween a physisorbed atom and a solid crystal surface is 
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determined by two main contributions. At large dis- 
tances the interaction is dominated by the attractive 
van der Waals polarization interaction, and close to the 
surface the interaction is repulsive, due in essence to 
the overlap of the electron wave functions of the two 
subsystems. At a distance z from a solid surface (the 
crystal occupying the lower half-space), a physisorbed 
atom is subject to a potential 

v{z)^vr{z) + vyw{z). (1) 

The van der Waals interaction is written as 

C 

Vvw{z) = --— -^/(fcc(z - 2:vw)), (2) 

\Z -2vw 7 

with the coefficient C given b y^^i^-*^ 

where a is the polarizability of the atom and e is the bulk 
dielectric function of the solid. The function / describes 
the fact that the van der Waals interaction saturates as 
z draws closer to the reference plane.— The vanishing 
of the response at short wavelengths introduces the cut- 
off wave vector /cc," for which we take the inverse of 
the Ps hard-core radius, = 1.9 We discuss Zvw, 

the so-called reference plane position, further on, when 
addressing the repulsive interaction. 

To obtain the dielectric function we require first to 
calculate the electronic structure of the crystal. We 
considered a-quartz (i.e., Si02 with trigonal structure, 
space group P3221), which is the most common and 
the one used in Ref. [ll|, and have employed the exper- 
imental structural parameters of Ref. |2^. We used the 
highly precise all-electron full-potential linearized aug- 
mented plane wave (FLAPW) implementation of density- 
functional theory.— The band contribution to the imagi- 
nary part of the dielectric function is calculated according 
to22 

e2{u:) = ^ hm VJ] l|(k + q,c|e^'i--|k,T;)p 

il q-tO ^-^ 
c,v k 

x5(ek+q,c - ek,« - (4) 

Here $7 is the crystal volume, and c and v denote the con- 
duction and valence bands, respectively. It is important 
to note that we use the screened-exchange local density 
approximation (sX-LDA) to the exchange-correlation po- 
tential. This approximation is known to give excellent 
results regarding the band gaps and optical properties 
of sp-semiconductors.^^ This is because of its superior 
description of the exchange-correlation hole-whose long 
range behavior is critical in semiconductors-compared to 
the well-known local density approximation (LDA) 2& In- 
deed, the sX-LDA band gap in Si02, -Eg = 8.8 eV, agrees 
very well with the measured value of 8.9 eV^^ while the 
LDA band gap is only of 6.1 eV. 
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FIG. 1: Imaginary part of the dielectric function of a-quartz 
for a polarization vector perpendicular to the c-axis. The con- 
tinuous line shows the total dielectric function, arising from 
interband and exciton excitations. To illustrate the excitonic 
effects, the dotted line shows the interband contribution only. 

Now, a recent ah initio study of excitonic effects in a- 
quartz has shown that these are crucial to understand 
the structure of its optical spectrum,— most notably the 
strong peaks observed near the absorption edge.'^^ We 
parametrize these contributions with a Hopfield term of 
the form ^jPj/{oJj — t^^ — iju;). In Fig. [T] we show 
both the band contribution and the total e2(w) for light 
polarized perpendicular to the c-axis of the crystal (a- 
quartz is optically uniaxial). Our result compares well 
with the aforementioned works. "^^ The real part of the 
dielectric function, obtained by the Kramers-Kronig re- 
lation, yields a calculated macroscopic dielectric constant 
Eoo = 2.38, which is in fact the experimental value. "^^ The 
dielectric function for imaginary argument is obtained by 
analytic continuation!^ Our result for the van der Waals 
coefficient is finally C = 13.87 eV ag. 

The polarizability of Ps is readily obtained by (re- 
duced) mass-rescaling that of hydrogen. For the latter, 
we take the parametrization of Dalgarno and Victor— 
a{uj) — J2n fn/i^n ^ '^ith the oscillator strengths 
and frequencies they give. 

The repulsive part of the interaction is essentially pro- 
portional to the valence electron density profile near the 
surface of the crystali^ We write analytically the repul- 
sive potential asi^ 

vn{z) = Voe-(^-^°)/', (5) 

where Vq is the intensity of the Ps work function, I is 
the electron density decay length, and zq is the so-called 
background edge position. For Vo we take the experimen- 
tal value, which is 1 eVM. To determine the two other 
quantities we calculate first the electronic density at the 
surface of the crystal. For this, we considered a single 
Si02 slab in which the crystal exposes the (Oil) face (the 
so-called AT-cut crystal) to the Ps atom, for direct com- 
parison with experiment. The slab thickness was ^ 24 
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1 2 3 4 curve) and of its sole excited state (dotted curve). The inset 

2 (^mj^^ shows the interaction potential. The potential saturates once 

it reaches the Ps work function value (1 eV). 



FIG. 2: Behaviour of the electron density at the surface of a 
a-quartz slab exposing a (Oil) surface, illustrated by a loga- 
rithmic plot of the electron density averaged over the section 
of the unit cell, as a function of the coordinate z perpendicular 
to the slab. 



ao, and the unit cell contained 6 formula units. Here, two 
technical points are worth mentioning. To eliminate spu- 
rious states and to avoid any possible dipole field effects 
on the top surface due to the finite thickness of the slab, 
the dangling bonds at the bottom surface of the slab 
were passivated with hydrogen at optimized positions. 
Further, the atoms near a surface generally reconstruct, 
changing their relative positions with respect to the in- 
finite crystal, and possibly affecting the electron density 
profile. Thus, the atom positions of the upper half of the 
slab were also relaxed and optimized. 

For given z, we define as the average of the valence 
electron density over a section of the unit cell parallel to 
the xOy plane. In Fig. Owe show a plot of Inn^, together 
with a linear fit. The origin of coordinates corresponds 
to the position of the uppermost layer of atoms (oxygen) . 
The behaviour becomes clearly linear at a certain point 
from the origin, defining both the electron density decay 
length, I, and the background edge position, zqM- The 
slope of the hnear fit yields = 2.06 a.u., and we find 
zq = 0.95 a.u. In our case, furthermore, zq also indicates 
the value of the reference plane position, z^w, for the 
van der Waals term. Indeed, in their work Zaremba and 
Kohn^^ found that in the case of insulators Zvw is located 
to a good approximation at the background edge. 

Note that the total potential v{z) saturates when 
reaching the work function value as z draws closer to 
the crystal. Thus, from this point on, we take a con- 
stant potential. Once the interaction potential is de- 
termined, we need to solve the Schrodinger equation 
[-{h'^/2m)\7'^ + v{z)]i}{z) = Eip{z) and see if it ad- 
mits any bound state solution. We solve this equation 



semi-analytically, exploiting an analytic expansion for 
z — > +00 of the form exp(— |_B|z) ^^^^p a„z^", a„ = 
— [Ca„-2 + n(n — l)a„_i]/2|i?|n. We find two bound 
states: a ground state ipo with energy Eq = —0.14 eV 
and a single excited state V'l with energy Ei — —5 meV. 
In Fig. [21 we show a plot of the probability density for 
both states, as well as of the potential v{z) in the inset. 
Thus, the Eq value we find is remarkably close to the 
—0.15 eV estimate by Sferlazzo and co-workers 

Given that there is no previous theoretical study along 
the above lines of the problem we are considering, and 
that there is little experimental information to date, 
it is useful to compare our results with some known 
facts regarding physisorption. First, we point out that 
there is an empirical, universal relationship between the 
van der Waals constant C, the macroscopic dielectric 
constant of the crystal, and the static atomic polariz- 
ability, put forward by Hoinkes,^'' which is expected to 
hold within a spread of roughly 30%. The relation is 
C = Ka{e-l)/{e+l), where K = 1.41 eV. With a = 4.5 
a.u. (static Ps polarizability) and e = 2.38 one finds 
C = 20.7 eV ttg, so that our value is 33% below. It is 
notable that the Hoinkes law still applies relatively well 
to such a light atom as Ps. Second, it is of interest to 
compare the energies we find with those of other systems. 
If we take H, which is the lightest normal atom, we find 
that its binding energy on graphite, for example, is 32 
meV (for other systems it is even lower). Thus, at first 
sight the ground state energy we find for Ps appears to 
be rather large. One can understand this, however, not- 
ing that the well depth in our case is ^ 1.7 eV, while in 
the case of H/graphite it is 43 meV. The fact that our 
value of Eq is such a small fraction of the well depth is 
of course due to the lightness of Ps. Also, the potential 
well in our case is so deep because of the strong polariz- 
ability of Ps (8 times that of H) and its relatively small 
work function. Thus, the values we find, although un- 
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usual, are overall well understood in the framework of 
physisorption. They appear to be robust also in the fol- 
lowing sense. For comparison, we have considered the 
case of the so-called Z-cut crystal, i.e., one exposing the 
(001) surface to the Ps atom. The results are very close 
to those of the AT-cut crystal, with a less than 1% dif- 
ference in the binding energies. Furthermore, the change 
in energies arising from surface reconstruction is of the 
same order of magnitude 

We now briefly discuss the possibility of Ps2 forma- 
tion, for which we invoke the Langmuir-Hinshelwood 
mechanismi^ Consider two Ps atoms initially ph- 
ysisorbed on the surface of a crystal. Since the Ps 
atoms attract each other with an effective Lennard- Jones 
potential,— they can eventually collide, recombine, and 
desorb. First, the fact that the atoms are trapped 
at the surface increases the chance that they do not 
simply rebound, and second, the Ps2 binding energy 
is — 0.44 eV. Thus, upon recombination the Ps2 
molecule releases internal energy, which is transfered to 
its z degree of freedom, and causes the molecule to be 
desorbed with an energy of ~ 0.16 eV (the initial total 
energy of the two trapped Ps atoms is 0.28 eV). 

In summary, we used a first principles approach to cal- 



culate the key parameters describing the Ps/a-quartz in- 
teraction potential at physisorption distances and showed 
that it indeed has a bound ground state, with an en- 
ergy very close to that estimated in experiment. This 
result, furthermore, is consistent with the possibility of 
Ps2 formation at the quartz surface via the Langmuir- 
Hinshelwood mechanism. We believe our results bring 
strong support to the idea of Ps trapping at the surface 
of certain insulators and to the fascinating prospect of 
creating many-Ps systems in this way. 
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